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The hybrid use of Loran-C with the Global Positioning 
System (GPS) has been shown capable of providing a sole- 
means of enroute air ladionavigabon. By allowing pilots 
10 fly direct to their destinations, use of this system is 
resulting in significant time savings and therefore fuel 
savings as well. However, a major error source limiting 
the accuracy of GPS is the intentional degradation of the 
GPS signal known as Selective Availability (SA). SA- 
induced position errors are highly correlated and far 
f x cff d all other error sources (horizcmtal position error! 
i(X) meters, 95%). Realtime mitigation of SA errors from 
the position solution is highly desirable. This paper 
discusses how that can be achieved. The stability of 
Loran-C signals is exploited to reduce SA errors. The 
theory behind this technique will be discussed and results 
using bench and flight data will be given. 


Introduction 

The hybrid use of Loran-C with the Global Positioning 
System (GPS) has been shown to be enable of providing 
a sole-means of enroute air radionavigation [1]. 

Standardixation committees such as the RTCA are 

currently working on developing minimum operational 
performance standards for this system. By allowing pilots 
to fly direct to their destinations, use of this system will 
result in significant time savings and therefore fuel savings 
as well. By not confining all aircraft to a small portion of 
the airspace (which results when using the Victor airways), 
the risk of collision undoubtedly will he reduced as well. 

However, a major error source limiting the accuracy of 
GPS is the intentional degradation of the GPS signal 
known as Selective Availability (SA). SA manifests itself 
in the form of erroneous orbital data broadcast by the 
satellites and in dithering of the satellite dock. The result 
is position determination which, according to the 
Department of Defense (DoD), will be in error by one 
hundred meters 95% of the time in the horizontal plane. 
Previous work performed at Ohio University showed that 
SA-induced position errors are highly correlated [2]. Since 
the correlation time is on the order of minutes, it follows 
that the error falls well within the passband of the 


Realtime mitigation of SA erron from the position 
solution is highly desirable. This paper discusses how that 
can be achieved. The stability of Loran-C signals is 
exploited to reduce SA errors. In the typical hybrid use of 
Loran-C and GPS, the Loian-C signal stability is not 
exploited. This stems from the reiativdy poor absolute 
accuracy of Loran-C (relative to GPS). However, it is 
possible to use the stability of Loian-C positioning to 
reduce SA-induced GPS positioning errors. The theory 
behind this technique will be discussed and results will be 
given. First, the phenomenon of SA will be described. 


Selective Availability 

As mentioned in the introduction, SA is an intentional 
corruption of the GPS signal by the DoD to limit the 
accuracy available to the public. The degradation is 
achieved in two ways. First, false satdlite orbit 
parameters are broad^t to the users. This results in 
incorrect positioning of the satellites in the navigation 
solution. SecotKlly, code and carrier tracking errors are 
induced through dithering the satellite clock (carrier 
frequency). The erroneous orbit paiaxneten lead to 
position errors which vary slowly thnwghout the satellite 
pass. Code-phase and carrier-phase errors due to the 
dithering of the satellite clock are random but also are 
highly ccnrelated. Correlation times of several minutes art 
typical. As a result, simple fllteiing schemes are not 
effective and aircraft will follow the deviations. Virtually 
all of the information available to date about SA has been 
gathered through data collection efforts by civilian 
organizations. Tlie DoD, however, has stated that SA 
shall be instituted in such a way as to yield horizontal 
position errors at a 95% level of ICX) meters [3], 


Mitigation MethodoLozY 

The heart of the mitigation scheme lies in the differences 
between Loian-C and SA-induced GPS position errors. 
Loian-C position erron in general are biased and noisy. 
The level of noise depends upon the receiver architecture 
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and spedficaliy upon the tracJong loop bandwidth. 
However, noise lev^ on the order of 5 to 10 meten can 
be achieved for airborne applications. The Loran-C 
posibon bias is primarily composed of unmodeled 
additional secondary phase factors (ASF). In general the 
bias does not remain constant over any given flight path 
but the variation is usually quite slow in comparison to the 
dock component of GPS SA error. This phenomenon is 
what mak» Loran-based SA mitigation possible. The 
long-term stability of the Loian-C measurements is 
exploited to smooth the SA-induced variations in the GPS 
measurements. 

Conceptually, the mitigation scheme works as follows. 
The Loran-C sensor computes the horizontal position of 
the aircraft. A vertical input is needed and is supplied by 
the barometric altimeter (again, a sensor which is biased 
but st^le). The combination provides a three- 
dimensional position of the aircraft. Range values are 
computed fiom the GPS satellites to the Loran- 
C/altimeter position. These range values are then the 
reference against which the measured GPS pseudoranges 
are filtered. 

Note that the technique depends upon the assumption that 
SA error is composed only of high frequency components 
relative to the Loran-C bias error variations. Strictly 
speaking, this assumption is not valid since the orbital 
component of SA error is slowly varying. However, as was 
shown in [2], the clock component of SA error has periods 
on the order of five to ten minutes. As such it is a high 
frequency error source relative to the non-noise 
component of Loran-C error. Although this has not been 
rigorously proven, flight data (to be shown later) supports 
the conclusion. Thus, the technique is able to reduce the 
dock component (or roughly speaking, the variance) of SA 
error. 

The filtering is accomplished by complementary Kalman 
filters which are applied to each pseudorange 
measurement [4,5] . The inputs to each filter are the given 
GPS pseudorange measurement and the corresponding 
range computed from the satellite to the Loran- 
dValtimeier pMtion. At each measurement epoch 
(current tune given by the index k), the complementary 
Kalman filter is executed as follows: 

- <-i ♦ (if Vj) 


P* - Pk-i * 9 (2) 
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where the subscript represents the time index. The 
supersci^ '-’and ’-t-’represent predicted and e-trimatwi 
quantities leqiectivdy. 'd*' represents the estimated 
pseudonnge with variance q. 'z’represenu the measured 
pseudorange with error variance r. Note that r is due 
primarily to SA. ’L'r^resenis the r^ge computed from 
the satellite to the I^rah-C/alSmaer position. 'p' 
r^tesents the prediction or estimation error variance, 'k' 
is the Kalman gain. In equation (1), the current 
pseudoiange prediction is computed by updating the 
previous pseudorange estimate with the difference 
between the current and previous Loran-C/altimeter 
ranges. The (nediction error variance is computed in 
equadon (2) and is used to compute the Kalman gain in 
equation Q). The difference between the measured and 
predicted pseudotanges is weighted by the Kalman gain in 
the computation of the current esdnuie (equadon 4). 
Finally, the current esdmadon error variance is computed 
(equadon S). 


Position Solution 

Given at least four GPS pseudoranges, position may be 
computed. As will be shown in the next two sections, 
significant reduction in SA-error may be achieved when 
using the mitigation technique just described. 

For both the simulation and flight test results (to be shown 
later), the ordinary least-squares (OLS) estimator is used 

to determin e p osition and clock bia s from the 

P*cudoran^. Li the absence of nv^imlwn'erTbn, the 
rela tion s h ip between satellite location, receiver location, 
clock bias and pseudorange is given by; 

- ^(x, -xf ♦(y, -y)* ♦ (z, -z)^ ♦ b (6) 

where R, is the pseudorange to the i* satellite, (Xj,yj,Zj)are 
the cowdinates of the satellite, (x,y,z)are the coordinates 
of the receiver and b is the receiver clock bias (converted 
to units of distance through multiplication by the speed of 
light). Since the receiver coordinates and clock bias must 
be solved for simultaneously, at least four measurements 
are required. 

However, instead of attempting simultaneous solution of 
non-iinear equations, the standard technique is to solve 
iteratively a set of equations which have been linearized 
about an initial estimated position and clock bias 
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(Xo,yo»Zo»bo)- This is achieved by forming a Taylor scries 
expansion and retaining the zeroth and first order terms: 
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The presence of measurement crron may be accounted 
for by the addition of an error vector: 

y = ^ c (ib) 

The ordinary least-squares solution is then given by: 

Ku " 


where is the range from the satellite to the initial 
position estimate. ^x, 5y, 6z and 6b represent the 
corrections to the initial estimates. If the initial estimate 
is close to the truth, no iterations are required. However, 
if the initial estimate is not dose, the corrections are used 
to update the initial estimate and the process is repeated. 
Convergence is declared if the magnitudes of the 
corrections are below a desired threshold. 


The partial derivatives are evaluated as follows: 
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Substitution of equations (8) through (11) into (7) yidds: 

an, » (ajt)au + (iy)«o * (5z)«u + (12) 


where: 


an, 


Rf - 


(13) 


Four pseudorange measurements allow for the following 
simultaneous set of equations: 
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which may be rewritten more succinctly: 
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(15) 


After one iteration then, the position and clock bias 
estimate is given by: 
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Simulatjon 

To determine the feasibility of the technique, a simulation 
was performed. A simple flight-path was modded with 
the aircraft traveling to the east for 900 seconds at 100 
meters/ second, followed by a 2g turn and then returning 
to the west (figure 1). For the sake of simplicity in the 
calculations, a static satdlite constdlation was modded. 
In order to focus on the effects of SA, all other GPS error 
sources were assumed to be zero. The Loran-C/altimeter 
errors were modeled in the position domain by a constant 
200 meter bias on each axis. 

The SA model was obtained from collected data using the 
System Identification procedure described in [2]. In order 
to modd SA rather than the combination of SA and 
receiver noise, integrated Doppler data (rather than 
pseudorange data) were used. The System Identification 
procedure yidded a 16th order autoregressive (AR) filter. 
When Gaussian white noise (of proper variance) is input 
to this filter modd, the output is statistically equivalent to 
the collected SA data. An example of the output is given 
in figure 2. 

The positioning errors resulting from the SA corruption 
arc given in figures 3 and 4. Both the east and north 
components of the position error exhibit similar 
characteristics to the SA error on the pseudorange 
measurements. As discussed earlier, the errors are highly 
correlated and reach up to 100 meters. However, use of 
the Loran-C/aldmeter data in the complementary Kalman 
filter yidds significant reduction of SA error (figures 5 and 
6). 
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Flight Test 

Although extremely encouraging » the simulation results 
were obtained using a simple model for Loran-C position 
errors. In order to verify Ac robustness of Ac technique, 
actual flight data was us^. This is necessary sAce Loran- 
C position error bias is spatially dependent. 

The flight data employed here were collected during a trip 
from Cleveland to AAcns, Ohio in Fall of 1990 (figure 7). 
It may be recalled that SA was temporarily turned off at 
Aat time because of military use of civilian GPS receivers 
during Operation Desert Shield [6]. As a result, Ae GPS 
horizontal posidonmg accuracy is on Ac order of 10-20 
meters [1]. For this flight, Ae GPS-derived position was 
Acrcforc used as a rough truA reference. 

SA was generated by Ac model described earlier and 
added to Ac raw GPS pscudorangc measurements (figure 
8). As expected, Ac Loran-C position error is biased but 
Ae bias is not constant wiA position (figures 9 and 10). 
As was done earlier, altimeter error was modeled as a 
constant 200 meter bias. Raw SA-mduced position errors 
arc as expected wiA large excursions and high correlation 
(figures 1 1 and 12). Agam, position errors after smoothing 
arc significantly reduced (figures 13 and 14). It is 
important to note Aat even in Ac face of spatially varying 
Loran-C position errors, Ae mitigation scheme continues 
to perform well. 


Conclusions 

A technique has been described whereby Ac stability of 
Loran-C signals arc exploited to reduce SA-induced GPS 
position errors. The viability of Ac technique has been 
confirmed using simulations as well as actual flight data. 
Future work will consider Ae possibility of realtime SA 
model identification. 
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Figure 8. Pscudorangc error due to SA. 
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Figure 9. Loran-C position error. 
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